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Abstract The bandstructure was calculated by the 
full-potential linearized augmented plane wave method. 
The result reveals two important insights to the novel 
second harmonic generation (SHG) of alpha-phase lithium 
iodate (a — L1IO3) crystal: the existence of finite intra- 
band momentum matrix elements due to the non-inversion 
symmetry of the crystal illuminating the potential of 
the intra-band transition, and the strong covalent bond- 
ing between the /-atoms and the ligand O-atoms re- 
sulting the condition of the double-resonance. An inter- 
band transition scenario in SHG as a — L1IO3 in nano- 
structure is proposed. The optical properties were cal- 
culated within the theoretical framework of the time- 
dependent perturbation of the independent-particle model. 
The dielectric tensors and the refractive index were eval- 
uated. Comparisons between the predictions and the re- 
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suits were made: the x-ray near edge absorption spec- 
tra; the refractive index at the static limit, and at finite 
frequencies. Possible factors attributing the calculation 
errors is discussed. 

PACS 42.70.Mp, 73.20.At, 78.20.Ci 



1 Introduction 

The alpha-phase lithium iodate (a — L1IO3), in a hexag- 
onal crystal as depicted in Fig. [1] 1 „2] , has been inten- 
sively studied in the past years due to its novel nonlinear 
optical [3], piezoelectric [HH], and acousto-optic [SBl] 
properties. Its nonlinear optical property, especially in 
the second harmonic generation (SHG), always receives 
the attentions from the electro-optical technology, which 
is recently rising up due to the all-optical information 
processing for modern telecommunications systems. Fur- 
thermore, very recently the development of a — LHO3- 
based nano-structural systems for the nonlinear optical 
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waveguide is very promising because of the preservation 
for the novel nonlinear optical functions in the bulk. [8]IS] 
In addition, the developments of photorefractive devices 
in real-time holography, and image processing, also ur- 
gently demand the ultraviolet photorefractive property 
of a — LHO3 to improve resolution and storage capacity. 
[TU] With the analogous characteristics: high refractive 
index, excellent optical transmittance in the visible and 
near-infrared region, and high dielectric constant, the 
newly invented room-temperature ferromagnetic semi- 
conductor, T1O2 in anatase structure under a Co-doping 
has inspired the spintronic technology toward the 
magneto-optic and the optoelectronic applications. Thus 
presumably, the a—LHO^ crystal should have the poten- 
tial to be developed as magnetically doped reformations 
or a associated substrate material in the same way. In 
spite of this material being invented very early, the the- 
oretical study at First-principles level on the electronic 
structure and the optical transition is surprisingly very 
limited. The resulting characters of bandstructure can 
further guide the searching and the tailoring of new non- 
linear optical materials. 

The present First-principles calculations were per- 
formed by the full-potential linearized augmented plane 
wave (FLAPW) method of density functional theory de- 
fined by the local density approximation (DFT-LDA). 
[I^[T^[T4lfT5lfT6l[T7] Testing by the f-sum rule QMS], the 
eigen-values and the momentum-matrix elements was 
verified to be valid to produce correct evaluations of op- 
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tical properties. By virtue of the analysis of the band- 
structure, the insights for the intra-band and inter-band 
transitions of the remarkable SHG property are discov- 
ered. The near-edge states dominating the optical tran- 
sition were further analyzed. The conclusion of it sup- 
plies an useful scenario in the classical dipole-oscillation 
framework to understand the dielectric response of the 
a — LHO3 crystal. In order to comparing with the exper- 
imental results about the bandstructure, a simulation of 
the x-ray near edge absorption spectra was incorporated 
here. 

On the calculation of optical functions, the time- 
dependent perturbation of independent particle model at 
the long-wavelength limit was applied, in which the eigen 
states of quasi-particle were approximated as those given 
by the DFT-LDA. In addition, in a non-perturbative way 
the time-dependent DFT , allowing for the ab-initio cal- 
culations on electron in external electromagnetic fields, 
also solves the analogous problem, so it is able to cap- 
ture the strong dynamic effect of system illuminated in 
the high power. [2"TU2"U] On the other hand, the recently 
well-developed Berry phase method [22] can also directly 
simulate the macroscopic polarization at the ab-initio 
level to the external field effect. Due to a greater de- 
mand of the computation, so far they are more spec- 
ified to the nonlinear optical simulations involved sig- 
nificant dynamic effect and the polarization in ferroelec- 
tric systems, respectively. However, the relation between 
the bandstructure information and the optical transi- 
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tion mechanisms, directly inspiring the tailoring the elec- 
tronic structure of materials, can be feasibly accessed by 
the perturbation method. Hence, plenty previous works 
applied it to calculate the dielectric and SHG tensors of 
transparent insulators. [^[^[^[^PTPSP^ISnilSTP^pg] 
So far comparing with experimental results, such a method 
makes good agreement for the dielectric response calcu- 
lations. 

The present article is outlined as follows. In Sec. [H 
the implementation of the FLAPW method and the for- 
malism for the dielectric response will be illustrated. In 
Sec. [H the resulting bandstructure and the dielectric 
functions will be exhibited. The insights for the novel 
SHG property and the optical transition structure will 
be discussed. In Sec.[U the comparisons with the experi- 
ments on the x-ray near edge absorption spectra, the re- 
fractive index at the static limit and the finite frequency 
will be given therein. Finally, a brief summary will be 
given in Sec. [5] 

2 METHODOLOGY AND FORMALISM 

2.1 FLAPW Method 

The modified FLAPW method, 'APW+lo' QS], was ap- 
plied via the implementation of the WIEN2K code. [TBI 
[T4l[T5l[T6] The exchange-correlation potential functional 
was defined by the generalized gradient approximation 
(GGA) parameterized by Perdew and Wang |34) . The 
core and the valence states were respectively calculated 
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relativistically and semi-relativistically. The muffin-tin 
radii were set to be 1.8, 1.83, 1.6 A for the Li, the J, 
and the 0-atom, respectively. The expansions of associ- 
ated Legendre polynomials for spherical harmonics of the 
wave function and of the non-spherical full-potential ex- 
pansion were truncated at I = 10 and I — 4, respectively. 
The parameter RK max was set to be 8.5. Additional lo- 
cal orbitals were added to incorporate low-lying valence 
states in the semi-core regime: 4d-states of the /-atoms, 
2s-states of the Li and of the O-atoms. The lattice con- 
stants a and c at energy minimum, 5.574 and 5.259, were 
obtained by the parabolic fitting of bulk-modulus calcu- 
lation, and the atomic coordinates: Li(0.0, 0.0, 0.07145); 
7(0.3333, 0.6666, 0.9991), and 0(0.2470, 0.3422, 0.8379) at 
the equilibrium were determined by the calculation of 
ionic relaxation. Those resulting lattice parameters for 
the present bandstructure calculation are very similar to 
the x-ray diffraction |T] reported previously. 

2.2 FORMALISM 

In the independent particle model, the optical conduc- 
tivity function a is expressed as |35j 

n,m 

5{ x )^^=-e^T)\ ( 2 ) 

where u n is the energy of the n-th band; Q denotes the 
volume of unit cell; u) is the photon energy, and p° TO is 
the a-component of the momentum matrix element in 
the Cartesian coordinate. Herein, the subscript n (to) 
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symbolizes a conduction (valance) band. The 6 func- 
tion was denned as a smooth Gaussian distribution with 
the r of 0.35 eV. Thus, the imaginary part of dielec- 
tric tensor Ei can be obtained by the relation: £%((*)) = 
(2tt /lo) a (w). By means of Kramers-Kronig transform, 
the real part e± and the refractive index n(uS) can be 
obtained as follows 

£1 m=i+8p r £ t" - dj , (3) 

Jo uj — w z 
„(«,) = ( V £ iH+^)+£lM)i/2, (4) 

The Brillouin-zone integration are achieved by the ir- 
reducible points of special-point sampling [36] . The con- 
vergence of the integration was tested. The error range 
was estimated to be at least less 0.1 percent by obtaining 
difference between two results given by the 12x12x11 
and the 23 x 23 x 21 mesh, respectively. 

3 RESULTS AND DISCUSSIONS 

3.1 BANDSTRUCTURE 

Firstly, the remarkable feature, all resulting bands within 
the range —6 to eV always appear in a number of pairs 
slightly split from one degenerate state, is distinguish- 
able in the Fig. [2l The resulting local density of state 
(LDOS), depicted in Fig. [3j indicates them to be de- 
rived by the O — O and I — O bonding. With respect to 
this point, according to the ' k-p' method, the band dis- 
persion about the -T-point can be approximated to the 
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second order in A k as [37 

E m ( k + A k) = E m ( k) + — A k- Pmm + h ^ ^ 

ft {A k ■ p mn )(A k ■ p nm ) 

+ m 2 ^ E m ( k)-E n ( k) ' U 

In a state with inversion symmetry, the parity of any its 
physical expectation is even, so the intra-band momen- 
tum matrix element p mm vanishes due to the odd parity 
of the momentum operator. Contrarily, in a state with- 
out the inversion symmetry, a finite p mm exists for the 
opposite situation. The former causes the first order term 
in the expansion to vanish so to give a perfect parabolic 
band curvature; however, the later results a non-zero lin- 
ear dependent energy-split added to the parabolic band. 
The scale of p mm should be much less than the fiA k, 
deduced from the resulting energy-split appearing to be 
very narrow. Then such a fine structure character can be 
viewed as the finger-print of the structural non- inversion 
symmetry of the a — LHO3 crystal. In fact, about the 
SHG the existence of finite p mm subjects the intra-band 
transitions [2"81l2T)] , though prohibited in the optical tran- 
sition, so the present work reveals the potential of this 
transition mechanism to generate certain contribution. 

Secondly, there is a large energy-split around lOeV 
resulted by the covalent bonding between the /-atoms 
and the ligand O-atoms according to the resulting LDOS, 
implies an extremely strong interaction associated with 
the bond forming. Actually it is even greater than the 
magnitude of the strong on-site Coulomb repulsion U in 
most of transition-metal or rare-earth atoms in the per- 



6 

ovskite crystals [38 | f39 ] , also widely applied as nonlinear 
optical materials 3 . Furthermore, the derived double- 
gap feature exhibited in the Fig. [2 separating the I — O 
bonding and anti-bonding states, as well as the inter- 
vening states localized on the O-atom near the Fermi- 
surface, naturally meets the isometric inter-band spac- 
ing condition for the double-resonance of the inter-band 
transition in SHG. However, the Li-atom was determined 
to make only little contribution to the aforementioned 
states. 

Predictably in the nano-structure system some sur- 
plus bands were induced within the gap because of the 
incorporation of surface localized states. According to 
the double-resonance scenario of the visual hole or the 
visual electron mechanism [53] illustrated in Fig. [JJ the 
original novel SHG in the bulk can be still preserved; 
however, those bands of surface localized states, unless 
they near the gap-edge, would be hard to give a signif- 
icant impact. On the other hand, after all presumably 
the number of the bands from surface localized states 
is much less than that from near-edge valence states in 
the bulk, also make it not be a major role in SHG. Such 
concepts should be useful to figure out the novel nonlin- 
ear optical properties in the bulk still preserving in the 
nano-structured systems. [U[S] 

The present calculation indicates the band gap as the 
type of allowed indirect transition. Based on the result- 
ing LDOS, the permitted transition between the gap- 
edge bands indeed occurs within each atomic muffin-tin 
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sphere, according to the angular momentum and parity 
selection rules of the atomic spectroscopy. The obtained 
indirect type agrees with the previous experimental con- 
clusion [JD] . It is deduced for the sake of the direction- 
ality of the near-edge p-states tending to maximize the 
band dispersion at the /"-point, and to minimize it at 
the zone corner. [IT] The obtained gap value, 3.8 eV, is 
slightly less than the measured edge-onset, ~ 4eV, in the 
absorption 42 . 43 , 40 , 44 and the transmission spectrum 

S3- 

The test of the f-sum rule QUITS], / °° ^-a(u)duj = 
^ifi = N e ff, was performed to exam the obtained band- 
structure and the momentum matrix elements. For the 
very less contribution, the (i-electrons of the /-atoms 
should be excluded, so the effective valence electrons 
is 52, being compatible with the results of N e tf: 56.83 
and 56.66 for the xx and the zz-component, respectively. 
This implies the obtained bandstructure and momentum 
matrix elements to be amenable to the following optical 
calculations. 

3.2 DIELECTRIC RESPONSE 

The resulting components of the dielectric tensor are ex- 
hibited in Fig. [5] In fact, the resulting yy-component is 
identity to the xx-cme, consisting with the experimental 
observations. The significant optical anisotropy behaves 
as the obviously different dispersions between the xx- 
and the zz-component. It is actually dominated by the 
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discrepancy of respective strength component, |p^ m | 2 , in 
the Eqn. ( [lj . Though the information of bandstrucure is 
hard to directly access the insight of this quantity, it 
still supplies the optical transition knowledge: the ini- 
tial state almost localizing on the O-atoms and the fi- 
nal state mainly derived from the I — O anti-bonding 
states. In fact, the absorption resonance-edge of the re- 
sulting imaginary part given by the gap separating the 
above two classes states, and the location of the absorp- 
tion peak, 7eV, equivalent to the energy-split between 
them, further solidate the aforementioned state informa- 
tion in the transition. Taking advantage of the classically 
electromagnetic radiation concept, the dynamical distri- 
bution of the charge-density in transition like a dipolc 
rapidly oscillating out of phase along the I — O bond, 
is useful to figure out the resulting optical anisotropy. 
It indeed gives an isotropic radiation on the xy-plane as 
the calculated result. 

4 COMPARISON WITH EXPERIMENTS 

The resulting X-ray absorption near edge spectra of the 
/-atom at the Lj, and the L/7/-edge are shown in Fig. [SI 
The used parameters of resolution identity to experimen- 
tal values: 0.75, and 0.66 eV |3S]; the inserted values of 
atomic natural widths are 3.46, and 3.08 eV [UJ for the 
Lj- and the L///-edge, respectively. The present calcu- 
lation agrees well with the previous measurements and 
their associated calculations of the multiple scattering 
theory [46], especially in the near edge regime. Such as 
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the experimentally observed white line feature at the 
Lj-edge, as well as the previously discovered pre-edge 
structure at the L/77-edge were reproduced. However, 
the discrepancy at the L/-edge might be due to the re- 
sulting 5j»-hybridization to be under-estimated. 

The resulting refractive indices and the correspond- 
ing experimental results are listed in Table [1] In fact, 
the adopted formula all are approximated at the long- 
wavelength limit, so the most adequate comparison to 
the experimental data, to drastically rule out the lo- 
cal field and the dynamic effects, should be right at the 
static limit. The constant, measured from a number of 
independent experiments [4"8"I|4"M5U] . is reproducible in 
the current calculation. Besides, the remarkable exper- 
imental negative birefringence of this crystal is also re- 
produced here. Since a slightly over-estimated ordinary 
and extra-ordinary component was produced here, this 
causes the disagreement with experimental values in the 
negative birefringence. Previous publications suggested 
the causes of discrepancy originating from the defect in 
details of the bandstructure [53"ll54| . so the unincorpo- 
rated non-local effect of the exact density functional [55] 
is deduced to make certain influence to them. 

The resulting dispersion of the ordinary n Q and the 
extra-ordinary refractive index n e at finite frequencies 
are shown in the Fig. [7j Generally, the present results 
perfectly match the frequency-dependent trend and give 
a consistent deviation with respect to the experimental 
measurements in the transparent regime. Some degree of 
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discrepancy from the experiments might be the result of 
un-incorporations of the aforementioned non-local effect 
and the dynamic many-body effect, and the local field 
effect. 

5 SUMMARY 

The truth of the bandstructure and the validity of the 
momentum matrix elements given by First-principles cal- 
culations are respectively verified by the comparison with 
the previous experimental result of X-ray absorption near 
edge spectra, and the test of the f-sum rule. The exis- 
tence of finite intra-band momentum matrix elements 
due to the non-inverse symmetry of the crystal is re- 
vealed by the resulting bandstructure, which dominates 
the intra-band transition in SHG. In addition, the energy- 
split of the I — O bonding results the condition of the 
double- resonance of the inter-band transition in SHG. 
The suggested scenario of the double-resonance in SHG 
as a — LHO3 in the nano-structure is an useful refer- 
ence to the tailoring of a same type of nonlinear optical 
systems. The present simulation, basing on the time- 
dependent perturbation in independent-particle model, 
can well capture the features of the linear dielectric re- 
sponse of the a — L1IO3 crystal. The discrepancies from 
the comparison with experiments are deduced from the 
ignored non-local effect in calculating the bandstructure, 
and the unincorporated many-body dynamic and local 
field effects in evaluating the dielectric functions. 
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Table 1 The refractive indices of a — LHO3 crystal at the 
static limit. The 'n \ the 'n e ' and the 'B' label the ordi- 
nary, the extra-ordinary component, and the birefringence, 
respectively. The A labels the operating wavelength. 
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Fig. 1 The unit cell structure of the a — LHO3 crystal. In 
the bonding sketch, each I-atom (purple color) locates at 
the top of pyramid based by three O-atoms (red color). The 
Li-atoms (grey color)reside the sites on the edge. 
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Fig. 2 The bandstructure of the a — LHO3 crystal. Here 
r : (0,0,0); M : 2ir/a(l/2, 0, 0); K : 2w/a(l/3, 1/3, 0), and 
A : 2tt/c(0, 0,1/2). 
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Fig. 3 The LDOS for the s- and p-states of the J-atom 
(top) and the O-atom (bottom). The inset of the top panel 
depicts the result of d-states of J-atom. 
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Fig. 4 The visual hole and the visual electron mechanims 
of the inter-band transition in SHG. The blue arrow depicts 
the original transition in the bulk, and the red specifies to 
the case for the bands of surface localized states only. 
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Fig. 5 The dispersion of the imaginary part, the real part, 
and the absolution of the dielectric tensor (in the unit of 
(eV ■ sec)" 1 ). The black and the red lines depict the e X x- and 
the £ zz -components of the dielectric function, respectively. 



Bandstructure and optical properties of a — L1IO3 crystal 
6 



13 



3 

N 



O 

2 



III 

L 1 -edge 


















#M&&>^ 1 1 1 







-20 -10 



10 20 30 40 50 
E(eV) 



3 



a 0.5 



o 




10 20 
E(eV) 



40 



Fig. 6 The X-ray absorption near edge spectra of the I- 
atom at the Li, and the Li u-edge. The present results, 
and the corresponding measurements; their associated simu- 
lations of the multiple scattering theory are respectively de- 
picted by the blue solid line, and the diamond symbols; the 
red dashed line. 
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Fig. 7 The dispersion of the refractive index. The solid line 
depicts the present result and the diamond, the cross, the 
square, and the triangle symbols exhibit the measurements 
of references [49], [51], [52], and [50], respectively. 



